Photon spectra from WIMP annihilation 
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If the present dark matter in the Universe annihilates into Standard Model particles, it must 
contribute to the fluxes of cosmic rays that are detected on the Earth, and in particular, to the 
observed gamma ray fluxes. The magnitude of such contribution depends on the particular dark 
matter candidate, but certain features of the produced photon spectra may be analyzed in a rather 
model-independent fashion. In this work we provide the complete photon spectra coming from 
WIMP annihilation into Standard Model particle-antiparticle pairs obtained by extensive Monte 
Carlo simulations. We present results for each individual annihilation channel and provide analytical 
fitting formulae for the different spectra for a wide range of WIMP masses. 

I. INTRODUCTION 



According to present observations of large scale structures, Cosmic Microwave Background (CMB) anisotropics and 
light nuclei abundances, the most important component of matter in the Universe cannot be accommodated within 
Oh the Standard Model (SM) of elementary particles. Indeed, Dark Matter (DM) cannot be made of any of the known 
particles, and this is one of the most appealing arguments for the existence of new physics. Indeed DM is a required 
' component not only on cosmological scales, but also for a satisfactory description of rotational speeds of galaxies, 
orbital velocities of galaxies in clusters, gravitational lensing of background objects by galaxy clusters, such as the 
Bullet Cluster, and the temperature distribution of hot gas in galaxies and clusters of galaxies. The experimental 
determination of the DM nature will require the interplay of collider experiments [l[ and astrophysical observations. 
These searches use to be classified in direct or indirect searches (see (2,01 f° r different alternatives). Nevertheless, 
non-gravitational evidence of its existence and a concrete understanding of its nature still remain elusive. Concerning 
direct searches, the elastic scattering of DM particles from nuclei should lead directly to observable nuclear recoil 
signatures. Although the number of DM particles which passes through the Earth each second is quite large, the weak 
interactions between DM and the standard matter makes DM direct detection extremely difficult. 
Q\ ■ On the other hand, DM might be detected indirectly, by observing their annihilation products into SM particles. 

Thus, even if WIMPs (Weakly Interacting Massive Particles) are stable, two of them may annihilate into ordinary 
matter such as quarks, leptons and gauge bosons. Their annihilation in different places (galactic halo, Sun, Earth, 
etc.) produce cosmic rays to be discriminated through distinctive signatures from the background. After WIMPs 
^ . annihilation a cascade process would occur. In the end the potentially observable stable particles would be neutrinos, 
gamma rays, positrons and antimatter (antiprotons, antihelium, antideuterions, etc.) that may be observed through 
different devices. Neutrinos and gamma rays have the advantage of maintaining their original direction thanks to 
their null electric charges. On the contrary, charged particles searches, such as those of positrons and other antimatter 
particles, are hindered by propagation trajectories. 

The detection of such indirect signals would not constitute a conclusive evidence for DM since the uncertainties 
in the specific DM interactions, DM densities and backgrounds from other sources are not fully understood yet. 
Nevertheless, this work precisely focuses on this kind of detection as an indirect method to get information about the 
DM nature, abundance and properties. 

Photon fluxes in specific DM models are usually obtained by software packages such as DarkSUSY and mi- 
crOMEGAs based on PHYTIA Monte Carlo event generator. In general, the total photon spectrum obtained from the 
addition of the contributions from different channels is obtained for the particular SUSY model under consideration 
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for a given WIMP mass. In this sense, it would be interesting to have a fitting function for the shape of the spectra 
corresponding to each individual annihilation channel and, in addition, determine the dependence of such spectra on 
the WIMP mass in a model independent way. This would allow to apply the results to alternative candidates for 
which software packages have not been developed, and obtain photon fluxes for arbitrary WIMP candidates. On the 
other hand, the information about channel contribution and mass dependence can be very useful in order to identify 
gamma-ray signals with specific WIMP candidates. 

The paper is organized as follows: in section II, we briefly review the standard procedure for the calculation of 
gamma-ray fluxes from WIMP pair annihilations. In section III, we comment on several aspects of detectors and 
backgrounds. Section IV is then devoted to the details of specific simulations performed with PYTHIA. In section 
V, we introduce the fitting formulae that will be used to describe the spectra and in section VI the results for the 
simulations, the fitted parameters and their dependence on the WIMP mass are presented. Then, in section VII 
we provide some information about the performed numerical codes obtained from our results and available online. 
Section VIII is then devoted to the main conclusions of the work. Finally, five appendices are provided in section IX 
to illustrate the obtained results for some studied annihilation channels. 



II. GAMMA RAY FLUX FROM DM ANNIHILATION 



Let us denote the DM mass by M and its thermal averaged annihilation cross-section into two SM particles (labelled 
by the subindex i) by (div) .Then the 7-ray flux from all possible annihilation channels is given by: 



Particle model dependent Dark matter density dependent 

where p is the DM density as a function of distance from its center r, which depends on the heliocentric distance s. 
The integral is performed along the line of sight (l.o.s.) to the target and averaged over the detector solid angle Af2. 

The first piece of the r.h.s. in ([1} depends on the particular particle physics model for DM annihilations. In 
particular, the self-annihilation cross sections is mainly described by the theory explaining the WIMP physics, whereas 
the number of photons produced in each decaying channel per energy interval involves decays and/or hadronization of 
unstable products, for instance quarks and gauge bosons. Consequently, the detailed study of these decay chains and 
non-perturbativc effects related to QCD is a hard task to be accomplished by any analytical approach. The second 
piece in ([TJ is a line-of-sight integration through the DM density distribution. We will discuss each of these pieces 
separately. 



d$ DM 



A. Particle Physics model 



Although annihilation cross sections are not known, they are restricted by collider constraints and direct detection. 
In addition, the thermal relic density in the range ficDM^ 2 = 0.1123 ± 0.0035 which is determined by fitting the stan- 
dard ACDM model to the WMAP7 data (Wilkinson Microwave Anisotropy Probe results for 7 years of observations) 
[H , the latest measurements from the B AO (Baryon Acoustic Oscillations) in the distribution of galaxies @ and the 
Hubble constant (i?o) measurement [6|, do not allow an arbitrary contribution from the DM gamma ray fluxes. 

As already mentioned, the annihilation of WIMPs is closely related to SM particle production. The time scale of an 
annihilation process is shorter than typical astrophysical scales. This fact implies that only stable or very long-lived 
particles survive to the WIMP annihilations and may therefore be observed by detectors. 

For most of the DM candidates, the production of mono-energetic photons is very suppressed. The main reason for 
such a suppression comes from the fact that DM is neutral. Thus, it is usually assumed that the gamma-ray signal 
comes fundamentally from secondary photons originated in the cascade of decays of gauge bosons and jets produced 
from WIMP annihilations. These annihilations would produce in the end a broad energy distribution of photons, 
which would be difficult to be distinguished from background. However, the directional dependence of the gamma ray 
intensity coming from these annihilations is mainly localized in point-like sources as will be discussed in the following 
section. This fact could therefore provide a distinctive signature. 

In conclusion, for a particular DM candidate, an unique annihilation channel may dominate, but in general, they 
all contribute. All those channels contributions produce a broad energy gamma ray flux, whose maximum constitutes 
a potential signature for its detection. Typically, this peak is centered at an energy that is one order of magnitude 
lower than the mass of the DM candidate. 



3 



On the other hand, a different strategy can be followed by taking into account the fact that the cosmic ray back- 
ground is suppressed at high energies. Primary photons coming from the Weicksacker- Williams radiation dominate 
the spectrum at energies close to the mass of the DM candidate and their signature is potentially observable as a 
cut-off ■ This approach has the advantage of being less sensitive to electroweak corrections which may be important 
if the mass of the DM candidate is larger than the electroweak scale [8J . 



B. DM density directionality 

The line of sight integration can be obtained from: 

(J) An = ^J An J ^ dn = ~mJ q m X d6 sin 9 1 ^ ds p2 (V s2 + a o _2sSoCOS0 ) (2) 

where 

J(V0= / dsp 2 (r). (3) 

./l.o.s. 

The angled brackets denote the averaging over the solid angle Af2, and s m i n and s max are the lower and upper limits 



of the line-of-sight integration: sq cos 6 ± y r\ — Sq sin 9. In this formula so is the heliocentric distance and r t is the 
tidal radius. 

Traditionally, the galactic center (GC) has attracted the attention of this type of directional analysis since standard 
cusped Navarro- Frenk- White (NFW) halos predict the existence of a very important amount of DM in that direction 
However, this assumption is in contradiction with a substantial body of astrophysical evidences [HI], and a 
core profile is not sensitive to standard DM candidates. On the contrary, cusped profiles are not excluded for the 
Local Group dwarf spheroidals (dSphs) that constitute interesting targets since they are much more dominated by 
DM. In this way, directional analysis towards Canis Major, Draco and Sagittarius or Segue 1 fl2j are more promising. 

An alternative strategy takes advantage of the large field of view of FERMI, that may be sensitive to the continuum 
photon flux coming from DM annihilation at moderate latitudes > 10°) [10]. Other proposed targets, as the Large 
Magellanic Cloud (LMC) [Ta |. are less interesting since their central parts are dominated by baryonic matter. 



III. DETECTORS AND BACKGROUNDS 



$max i n Eq. @ is the angle over which we average, and is bounded from below by the experimental resolution of 
the particular detector: 

An = 27!-/ d0sm0 = 27r(l-cos(0 inax )). (4) 
Jo 

The quoted point spread function widths for the various experiments are typically: 0.4° (EGRET), 0.1° (CANGAROO- 
III, FERMI, HESS, MAGIC and VERITAS). EGRET and FERMI are satellite detectors with low energy thresholds 
( about 100 MeV), high energy resolution (~ 15%) but only moderate angular precision. The others are atmospheric 
Cerenkov telescopes (ACTs) with higher thresholds (w 100 GeV) but better angular resolution. Typical reference 
sizes for the solid angle are AO = 10~ 5 sr for ACTs and FERMI and Ail = 10~ 3 sr for EGRET. 

There are different main sources of background for the signal under consideration: hadronic, cosmic-ray electrons, 
localized astrophysical sources and the diffuse 7-rays. The latter is negligible for ACTs, but only the last two are 
present for satellite experiments like FERMI or EGRET. 

For heavy WIMPs, the produced high-energy gamma photons could be in the range 30 GeV-10 TeV, detectable 
by ACTs such as HESS, VERITAS or MAGIC. On the contrary, for lighter WIMPs, the photon fluxes would be in 
the range detectable by space-based gamma ray observatories [lj] such as EGRET, FERMI or AMS, with better 
sensitivities around 30 MeV-300 GeV. 



IV. MONTE CARLO SPECTRA GENERATION: TECHNICALITIES 



In this section, we explicitly specify how gamma ray s spectra have been generated. We have used a widely known 
particle physics software, PYTHIA (version 6.418) [15|], to obtain the results we are about to present. In a first ap- 
proximation, the WIMP annihilation is described by two separated processes: The first one describes the annihilation 
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Channel \ Mass (GeV) 


100 


125 


150 


200 


250 


350 


500 


1000 


W+W~ 


5.21 




1.91 


6.85 




7.83 


2.91 


2.85 


zz 


0.42 


6.01 


2.91 


14.9 




14.2 


2.81 


2.02 


ti 








0.70 


0.86 


0.32 


2.81 


1.41 



Table I: Total number of photons - in 10 7 units — generated from W + W , ZZ and tt channels for different WIMP masses. 
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200 
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1000 
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5- 10 4 


r + r - 


2.25 


2.25 


2.23 


1.07 


2.81 


2.33 


8.41 


7.80 




2.25 


2.25 


2.23 


1.07 


2.81 


2.33 


8.41 


7.80 



Table II: Total number of photons - in 10 7 units - generated from t+t and u channels for different WIMP masses. 

of WIMP particles and its output which are particle-antiparticle SM pairs. The details are contained in the theory 
describing the WIMP physics. The second process considers the evolution (decays and/or hadronization) of the SM 
unstable products, for instance, quarks and gauge bosons. Unfortunately, a first-principle description of this latter 
step is too complex due to chain decays and non-perturbative QCD effects. 

As we mentioned above, in this work we have used PYTHIA to generate the photon energy spectra starting from 
pairs of SM particles, where each pair respects WIMP annihilation quantum numbers like neutral charge and color 
singlet. As will be described below, we will allow for final state radiation from charged particles to contribute to 
the photon spectra. Due to the expected velocity dispersion of DM, we expect most of the annihilations to happen 
quasi-statically. This fact offers the center of mass (CM) frame as the most suitable frame to produce the photon 
spectra. Hence, the process is described by the total energy: 

Ecu — 2 M (5) 

where M is the mass of the WIMP particle. Therefore, by considering different CM energies for the SM particles 
pairs in each WIMP annihilation process we are indeed studying different WIMP masses. The procedure to obtain 
the photon spectra is thus straightforward, except for the particular case of the t quark. For any given pair of SM 
particles which are produced in the WIMP annihilation, we count the number of photons in each bin of energy and 
then normalize them to the total number of simulated pair collisions. The bins which we have considered in the x 
variable, x = E^/M, are: [1CT 5 , 10" 3 ], [1(T 3 , 0.2], [0.2, 0.5], [0.5, 0.8] and [0.8, 1.0]. Nevertheless, for some studied 
channels more precision was needed in some particular energy intervals and additional bins were considered. 

The number of simulated collisions in each bin was fixed a priori but it was changed if required in order to provide 
suitable statistics in the number of produced photons. For instance, for the high energy bins many collisions are 
required to get a significant number of photons, whereas for low-intermediate energy, many photons are usually 
produced even for a small number of collisions. The total number of photons corresponding to the different generated 
pairs in terms of the WIMP mass are presented in Tables UJ [TT] and Mil These results will be plotted in the Appendix 
E, Figure 9 at the end of the paper. 

The SM particle pairs decays generated are W and Z gauge bosons, r and fi leptons and u, d, s, c, b and t quarks. 
For each annihilation channel we have studied the gamma ray spectra produced for different WIMP masses. The 
result of the simulations were fitted to analytical expressions as is described in the following section. 



A. Final state radiation 

If the final state in the annihilation process contains charged particles, there is a finite probability of emission of an 
additional photon. This is discussed in detail in [l6| . In principle there are two types of contributions: that coming 



Channel \ Mass (GeV) 


50 


100 


200 


500 


1000 


2000 


5000 


7000 


8000 


uu 


2.05 


11.9 


2.42 


2.81 


3.82 


10.8 


5.91 




2.11 


dd 


1.04 


1.96 


2.42 


2.81 


2.81 


2.81 


2.31 






ss 


15.3 


2.00 


1.97 


2.81 


9.82 


2.71 


2.71 


11.0 




cc 


2.41 


1.99 


16.8 


2.81 


2.81 


3.81 


12.0 




3.00 


bb 


11.7 


1.91 


2.62 


2.61 


8.81 


2.20 


3.81 




1.70 



Table III: Total number of photons - in 10 7 units - generated from uu, dd, ss, cc and bb channels for different WIMP masses. 
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from photons directly radiated from the external legs, which is the final state radiation we have considered in the work, 
and that coming from virtual particles exchanged in the WIMP annihilation process. The first kind of contribution can 
be described for relativistic final states by means of an universal Weizsacker- Williams term fundamentally independent 
from the particle physics model 16]. On the other hand, radiation from virtual particles only takes place in certain 
DM models and is only relevant in particular cases, for instance, when the virtual particle mass is almost degenerate 
with the WIMP mass. Even in these cases, it has been shown (T3] that although this effect has to be included for 
the complete evaluation of fluxes of high energy photons from WIMP annihilation, its contribution is relevant only 
in models and at energies where the lines contribution is dominant over the secondary photons. For those reasons 
and since the aim of the present work is to provide model independent results for photon spectra, only final state 
radiation was included in our simulations. 



B. The case for t quark decay 



The decay of top quark is not explicitly included in PYTHIA package. We have approximated this process 
by its dominant SM decay, i.e. each (anti) top decays into and (anti) bottom. In order to maintain any 

non-perturbative effect, we work on an initial four-particle state composed by W + b coming from the top and W~b 
from antitop, which keeps all kinematics and color properties from the original pair. Starting from this configuration, 
we have forced decays and hadronization processes to evolve as PYTHIA does and therefore, the gamma rays spectra 
corresponding to this channel have also been included in our analysis. 



V. ANALYTICAL FITS TO PYTHIA SIMULATION SPECTRA 



In this section we present the fitting functions used for the different channels. According to the PYTHIA simulations 
described in the previous section, three different paramctrizations were required in order to fit all available data from 
the studied channels. The first one for quarks (except the top) and leptons. Then, a second one for gauge bosons W 
and Z and a third one for the top. 



A. Quarks and leptons 



For quarks (except the top), r and (i leptons, the most general formula needed to reproduce the behaviour of the 
differential number of photons per photon energy may be written as: 



x L - b — — = aiexp -b 2 x " 2 -5- + — =- +p lnpl-i (6) 

ax V x dl x d2 1 x 

In this formula, the logarithmic term takes into account the final state radiation through the Weizsacker- Williams 
expression [HI Il8j|. Nevertheless, initial radiation is removed from our Monte Carlo simulations in order to avoid 
wrongly counting their possible contributions. 

Strictly speaking, the p parameter in the Weizsacker- Williams term in the previous formula is (M / ' m part icie) 2 where 
^particle is the mass of the charged particle that emits radiation. However in our case, it will be a free parameter 
to be fitted since the radiation comes from many possible charged particles, which are produced along the decay 
and hadronization processes. Therefore we are encapsulating all the bremsstrahlung effects in a single Weizsacker- 
Williams-like term. 

Concerning the y, lepton, the expression above (j6|) becomes simpler since the exponential contribution is absent. 
The [i~ decays in e~D e v^ with a branching ratio of ~ 1 and therefore the only contribution in addition to its own 
bremsstrahlung, is provided by the radiation coming from the electron. The total gamma rays flux is thus well fitted 
by: 

x i.s^2. = qx ^\n\p(l-x 1 )} x2 ~ 2x + 2 (7 ) 
ax x 

where the I parameter in the logarithm is needed in order to fit the simulations as will be seen in the corresponding 
sections. 
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Let us mention at this stage that for the gamma rays obtained from electron-positron pairs, the only contribution 
is that coming from bremsstrahlung. Therefore, the previous expression ([7]) is also valid with q — «QED/' r ! P — 
(M/m e -) 2 and I = 1. This choice of the parameters corresponds of course to the well-known Weizsacker- Williams 
formula. 



B. W and Z bosons 



For the W and Z gauge bosons, the parametrization used to fit the Monte Carlo simulation is: 

This expression differs from the expression ([5]) in the absence of the additive logarithmic contribution. Nonetheless, 
this contribution acquires a multiplicative behaviour. The exponential contribution is also quite simplified with only 
one positive and one negative power laws. Moreover, ai, n\ and q parameters appear to be independent of the WIMP 
mass M as will be seen in the corresponding section. The rest of parameters, i.e., b±, ci, di, p and j, are WIMP 
mass dependent and will be determined for each WIMP mass and for the W and Z separately. In both cases we have 
covered a WIMP mass range from 100 to 10 4 GeV. Nonetheless, at masses higher than 1000 GeV, we have observed 
no significant change in the photon spectra for both particles. 



C. t quark 



Finally, for the top, the required parametrization turned out to be: 



Likewise the previous case for W and Z bosons, gamma-ray spectra parametrization for the top is quite different 
from that given by expression ([6]). This time, the exponential contribution is more complicated than the one in 
expression ([5]), with one positive and two negative power laws. Again, the additive logarithmic contribution is absent 
but it acquires a multiplicative behaviour. Notice the exponent I in the logarithmic argument, which is required to 
provide correct fits for this particle. 

The covered WIMP mass range for the top case was from 200 to 10 5 GeV. Nevertheless, at masses higher than 
1000 GeV we have observed again that there is no significant change in the gamma-ray spectra. 



VI. RESULTS FROM PYTHIA SIMULATION 



In this section we present the results of our fit of the parameters given by expressions (HJ), ([8]) and ([9]) after 
having performed the PYTHIA simulations described in section IV. For each studied channel, we have considered the 
possibility of parameters depending on the WIMP mass. 

Once the parameters in expressions ©, © and © have been determined for each channel and different WIMP 
masses, it is possible to study their evolution with the WIMP mass M. Some parameters in expressions are WIMP 
mass independent and take values that depend on the studied channel. The rest are WIMP mass dependent. 

For some channels and in some range of WIMP masses, we observed that this dependence was given by a simple 
power law. In fact, for a given channel (i) and a generic mass dependent V parameter, a simple power-law scaling 
behavior would correspond to an expression like 

pW(M) = m pl ,)MM-) (10) 

with m V (i) and n-pm constant values to be determined for the different studied channels. Values of m V (i> and 
and their range of validity are presented for each studied channel in the following. 
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WIMP mass (GeV) 


bi 


Cl 


di 


V 


j 


100 


9.48 


0.651 


0.292 


973 


0.790 


-LiJU 


8 87 


U.OUO 






n qi q 


200 


8.64 


0.882 


0.250 


684 


0.955 


350 


8.56 


0.907 


0.245 


593 


0.991 


500 


8.51 


0.917 


0.244 


560 


0.996 


1000 


8.45 


0.931 


0.242 


535 


1.000 



ai = 25.8 ; m = 0.510 ; q = 3.00 



Table IV: W boson: b\, ci, di, p and j parameters corresponding to JHJ in the W + W channel for different WIMP masses. Mass 
independent parameters in ((Sj for this channel are presented at the bottom of the table. 



Parameter 


WIMP mass interval (GeV) 


Fitting power law(s) 


bi 


100 < M < 200 
200 < M < 1000 


0.0433M u ' bb + 46.4M _u - 3t " 
9.29M- 00139 


Cl 


100 < M < 200 
200 < M < 1000 


-27M- U -' ZW + 35.0A_T uub43 - 16.5 
0.743M 0331 


di 


100 < M < 240 
240 < M < 1000 


2.64 ■ 10 _4 M 1 - U3 + 2.28M- U ' 4YU 
0.265M- 00137 


p 


200 < M < 1000 


lO^M" 1 - 13 +285M UU " J4 


j 


385 < M < 1000 


0943M u.(ju« M 



Table V: Parameters corresponding to (1 1 P for W boson. It can be seen that p parameter follows two different power laws depending on 
the WIMP mass interval. For the remaining mass dependent parameters there is a unique power law behavior in the WIMP mass interval 
350 < M < 1000. 



A. W boson 

As commented above, the correct parametrization for the W boson simulations was given by expression (|8]). For 
this boson, there are five mass-dependent parameters: bi, ci, d\, p and j whose values are detailed in Table llVl The 
mass independent parameters are a\ = 25.8, n\ — 0.51 and q — 3.00. The mass range considered for this boson is 100 
to 10 5 GeV. In fact, from M — 1000 GeV, the photon spectrum does not change. The parameters obtained fit the 
enegy spectra from x = 2 • 10 -4 till the end of the allowed interval. It can be seen that for low masses the spectrum 
does not end at x = 1 but at smaller energies (e.g. x ~ 0.78 for M = 100 GeV) and as masses get higher, the energy 
tail approaches x = 1. 

Some of these results are presented in Figure 1 in Appendix A for four WIMP masses: 100, 200, 350 and 1000 
GeV. Besides, mass dependent parameters b\, c\, d%, p and j were presented in the same Appendix in Figure 2. 

Concerning the scaling behavior of these mass dependent parameters given by expression (|10[) . we obtain that bi, 
Ci and j parameters scale with a simple power law of M at high masses. In fact, b\ and c\ parameters follow a two 
power-law behavior at low masses. For d\ parameter, we find that the sum of two power laws covers this high masses 
interval, whereas a simple power law at low masses is obeyed. Parameter p scales with two power laws in the whole 
studied mass interval. These results are shown in Table [V] 



B. Z boson 

For the Z boson the correct parametrization is again the one given by expression (|8]). For this boson there are five 
mass-dependent parameters: bi, c\, d\,p and j which are detailed in Table IVT1 The mass independent parameters are 
ai = 25.8, n\ ~ 0.5 and q = 3.87. The studied WIMP mass range for this boson was from 100 to 10 5 GeV. However, 
above M = 1000 GeV the energy spectrum does not change as can be seen from our simulations. 

The chosen parameters values fit the photon spectra from x = 5 ■ 10 -4 till the end of the allowed interval. As for the 
W case, it can be seen that for low masses the spectrum does not end at x = 1 but at smaller energies (e.g. x ~ 0.7 
for M = 100 GeV) and as masses get higher, the high-energy tail approaches x = 1. 

Concerning the power-law scaling of the parameters with M, we obtained that parameters b\, c\, d\ ad j follow 
a simple power-law behavior for high WIMP masses. Parameter p follow a two sum power-law behavior for masses 
higher than 170 GeV. Concerning d\ parameter, the whole accessible WIMP mass interval is covered by different 
either one or two power laws. These results can be seen in Table IVlTl 
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WIMP mass (GeV) 


bi 


Cl 


rfi 


P 


3 


100 


10.3 


0.498 


0.323 


7010 


0.702 


125 


9.74 


0.612 


0.294 


4220 


0.836 


150 


9.49 


0.675 


0.280 


3850 


0.894 


200 


9.28 


0.734 


0.268 


3630 


0.943 


350 


9.02 


0.800 


0.257 


3380 


0.978 


500 


8.95 


0.813 


0.255 


3260 


0.988 


1000 


8.91 


0.819 


0.254 


3140 


0.997 



ai = 25.8 ; m = 0.5 ; g = 3.87 



Table VI: Z boson: 6i, ci, di, p and j parameters corresponding to (|8j in the ZZ channel for different WIMP masses. Mass independent 
parameters in $5$ for this channel are presented at the bottom of the table. 



Parameter 


WIMP mass interval (GeV) 


Fitting power law(s) 


6i 


500 <M< 1000 


9 36M -u.uuviu 


Cl 


465 < M < 1000 


0.765M ,U,UUB0 


di 


100 < M < 191 
191 < M < 360 
360 << M < 1000 


0. 00999 A/ UMU +21.5M" 1 '" 1 
2.02 • 10- 9 M 2 ' 56 + 0.491AT 115 
0.272M- 00990 


P 


170 < M < 1000 


8550A/- Uibb + 0.476Af uy84 


3 


350 < M < 1000 


0.884A/ u ' ulYb 



Table VII: Parameters corresponding to (1 1 I t for the Z boson. It can be seen that all mass dependent parameters for the Z boson follow 
a simple power-law scaling at intermediate and high masses. 



C. t quark 

For the top, there are six mass dependent parameters: b%, m, C2, p q and I which are detailed in Table [Villi The 
mass independent parameters are a± = 290, c\ = 1.61, d\ = 0.19 and c?2 = 0.845. The mass range for this quark is 
from 200 to 10 5 GeV. Nevertheless, from 1000 GeV onwards, the photon spectra do not change as was proven by 
considering several higher masses. The chosen parameters fit the spectra from x = 10 till the end of the allowed 
interval. Again for low masses, the spectra do not end at x = 1 but at smaller energies (e.g. x ~ 0.7 for m — 200 GeV 
) and, as masses get higher, the spectral tail approaches x = 1. 

Some of these results are presented graphically in Figure 3, Appendix B for four WIMP masses: 200, 250, 500 and 
1000 GeV . Also in this Appendix, mass dependent parameters bi, n\, C2, p, q and I are plotted in Figure 4. 

Concerning the scaling behavior of the C2, p, q and I parameters, they obey a simple power law in the whole 
accessible WIMP mass range. Nevertheless, for b\ and Ci parameters the simple power law behavior starts from 
masses bigger than 350 GeV. These results can be seen in Table IIXI 



WIMP mass (GeV) 


bi 


rii 


C2 


P 


<7 


I 


200 


14.4 


0.477 


3.34 ■ 10" 4 


1.34 


1.76 


4.42 


250 


13.5 


0.457 


1.54 ■ 10~ 4 


1.95 


1.96 


4.14 


350 


13.0 


0.448 


5.99 • nn 6 


3.78 


2.32 


3.74 


500 


12.8 


0.442 


1.69 ■ 10~ 5 


7.40 


2.75 


3.36 


1000 


12.4 


0.436 


1.80 ■ 10" 6 


30.0 


3.85 


2.72 



ai = 290 ; ci = 1.61 ; di = 0.19 ; d 2 = 0.845 



Table VIII: t quark: bi, m, C2, p, q and I parameters corresponding to 1(9} in the tt channel for different WIMP masses. Mass independent 
parameters in (O for this channel are presented at the bottom of the table. 
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Parameter 


WIMP mass interval (GeV) 


Fitting power law(s) 


bi 


200 < M < 350 
350 < M < 1000 


9.32M + 11.0 ■ 10 M 

i a a n f— 0.0400 

16.4M 


m 


200 < M < 300 

oUU \ IVl \ 1UUU 


21AM + 0.00867A1 

n KKQA/f- ' 0379 

u.ooyiVi 


C2 


200 < M < 1000 


sgioAf"' 5 ^ 


p 


200 < M < 1000 


5.78 ■ 10~ b M i SH 


Q 


200 < M < 1000 


0.133i\/ U4S8 


I 


200 < M < 1000 


21.9A/" U " ,U ^ 



Table IX: Parameters corresponding to (1 1 P for the t quark. It can be seen that all mass dependent parameters for the t quark follow a 
simple power-law scaling behavior at intermediate and high WIMP masses. Parameters b\ and ri\ follow the simple power law for M > 300 
GeV whereas the rest of parameters presented do so in the interval 200 < M < 1000. 



WIMP mass (GeV) 


rti 


P 


25 


10.1 


221 


50 


10.0 


767 


100 


9.91 


2520 


200 


9.80 


8660 


500 


9.67 


4.01 • 10 4 


1000 


9.57 


1.35 • 10 5 


10 4 


9.25 


4.80 • 10 6 


5- 10 4 


9.14 


5.44 • 10 7 



oi = 14.7 ; 6i = 5.40 ; b 2 = 5.31 ; n 2 = 1.40 ; 
ci = 2.54 ; di = 0.295 ; c 2 = 0.373 ; d 2 = 0.470 ; q = 0.00260 



Table X: r lepton: n\ and p parameters corresponding to J5J in the r + r channel for different WIMP masses. Mass independent 
parameters in (|6j for this channel are presented at the bottom of the table. 

D. Leptons and quarks 

For the rest of the quarks and leptons, the parametrization given in (j6|) is completely valid. Now we present results 
for t and fi leptons and all quarks except for the top. 

1. t lepton 

For the r lepton, there are only two mass dependent parameters in the spectra fitting function ©: n\ and p. The 
remaining parameters are mass independent for this particle and their values are a\ = 14.7, b\ — 5.40, hi — 5.31, 
n 2 = 1.40, ci = 2.54, d x = 0.295, c 2 = 0.373, d 2 = 0.470 and q = 0.00260. These results are presented in TableKl 

In this case, the WIMP mass interval considered ranges from 25 to 5 • 10 4 GeV. For masses higher than 5 • 10 4 GeV, 
the spectra do not seem to change, within the statistical uncertainties, with respect to that corresponding to 5 • 10 4 
GeV. 

The ni parameter scales with the WIMP mass as a simple power law for M < 5 • 10 4 GeV. For the other mass 
dependent p parameter, the power-law behavior is valid in two separated intervals with an inflection point in the 
behavior at M = 1000 GeV. These results can be seen in Table IXll 

Some of these results are presented graphically in Figure 5, Appendix C for four WIMP masses: 25, 100, 1000 and 
5 • 10 4 GeV. Also in this Appendix, mass dependent parameters rii and p are presented in Figure 6. 

For this particle, it is worth mentioning the increasing contribution of the logarithmic term in (j6]) as the WIMP 
mass increases. This fact can be seen in the presented plots from x = 0.5 onwards. As a consequence, the values of p 
parameter increase as the WIMP masses increase. 



2. jj, lepton 

For the fi particle and according to expression ([7]), there are only three mass dependent parameters: q, p and 
These values are presented in Table IXII1 In this case, the considered range for WIMP masses is from 25 to 5 • 10 4 
GeV. 
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Parameter 


WIMr mass interval (GeVJ 


Fitting power law(s) 


m 


25 < M < 10* 
10 4 < M < 5 ■ 10 4 


10.6M- - 0148 
-7.00M" 1 ' 99 + 179M" ' 763 + 9.09 


P 


25 < M < 1000 
1000 < M < 5 • 10 4 


0.773M iYb 
3 07M i.55 



Table XI: Parameters for expression l|10|l for r lepton. It can be seen that ni and p parameters follow power-law behaviors. 



WIMP mass (GeV) 


P 


q 


I 


25 


9510 


3.37 


10" 


■i 


0.787 


50 


23600 


3.40 


10" 


-3 


0.642 


100 


54600 


3.45 


10" 


3 


0.579 


200 


1.12 • 10 5 


3.50 


10" 


-3 


0.548 


500 


2.54 • 10 5 


3.61 


10" 


3 


0.523 


1000 


4.13 ■ 10 5 


3.70 


10" 


-3 


0.511 


5000 


1.18 • 10 6 


3.91 


10" 


-3 


0.484 


10 4 


1.84 • 10 6 


4.00 


10" 


3 


0.474 


5- 10 4 


5.29 • 10 6 


4.23 


10" 


-3 


0.454 



Table XII: fi lepton: Parameters corresponding to (0 for different WIMP masses. All parameters in expression (0 are WIMP mass 
dependent. 

The scaling of the p parameter with the WIMP mass shows two well differentiated regimes, with different asymptotic 
power laws: one from M = 25 GeV to M = 100 GeV, and another from M = 750 GeV to M = 5 • 10 4 GeV. On the 
other hand, q and / parameters present a sum of two power laws evolution in the whole studied WIMP mass range. 

As for the r lepton, the flux of photons increases as the WIMP mass increases. In this case, the q parameter 
increases as the WIMP masses do so, instead of the p parameter as was the case for the r. 



3. u quark 

The mass independent parameters are a\ = 5.58, 62 = 5.50, c\ = 0.315, C2 = 0.0 (therefore di is irrelevant) and 
q = 9.30 • 10~ 4 . The mass dependent parameters are 61, m, n,2, d\ and p. These results are presented in Table IxTVl 
The analyzed mass range for this quark is from 50 to 8000 GeV. 

The spectra of the two highest studied masses (5000 and 8000 GeV) clearly differ in the low energy interval. 
Therefore no conclusion can be made about the existence of an asymptotic high masses limit in the spectral shapes. 
Concerning the mass evolution of the parameters for this quark, we observe simple power-law behaviors for both b\ 
and p parameters in the whole studied WIMP mass range interval. On the other hand, tlx, ni and d\ parameters are 
fitted by a sum of two power laws in the studied range. These resuls can be seen in Table IXVl The chosen values for 
the parameters turn out to fit the spectra from x = 5 • 10~ 4 till the end of the allowed energy interval. Nevertheless, 
for some masses, the fit also applies for lower energies, i.e. lower x values, up to 10~ 4 . 



Parameter 


WIMP mass interval (GeV) 


Fitting power law(s) 


P 


25 < M < 100 


176M 1 -^ 




750 < M < 5 ■ 10 4 


4530Af 653 


q 


25 < M < 5 • 10 4 


0.00230A^r u - yil + 0.00291M u - ua4S 


I 


25 < M < 5 • 10 4 


0.626M" uu ' iuu + 16.4M- 1 " 44 



Table XIII: Parameters for expression l|10| l for fi lepton. It can be seen that p parameter follows simple power laws in two different WIMP 
mass regimes. Nevertheless, for q and I parameters a sum of two power laws accounts for the WIMP mass dependence in the whole studied 
WIMP mass interval. 
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WIMP mass (GeV) 


bi 


m 


n 2 


di 


P 


50 


3.60 


2.77 


0.585 


0.383 


129 


100 


3.75 


2.64 


0.551 


0.355 


225 


200 


3.88 


2.54 


0.521 


0.332 


409 


500 


4.04 


2.44 


0.490 


0.308 


856 


1000 


4.18 


2.40 


0.472 


0.293 


1540 


2000 


4.34 


2.39 


0.463 


0.281 


2800 


5000 


4.55 


2.38 


0.450 


0.266 


6000 


8000 


4.67 


2.34 


0.448 


0.259 


8900 



|ai = 5.58 ; b 2 = 5.50 ; a = 0.315 ; c 2 = 0.0 ; q = 9.30 ■ 10~ 4 .| 

Table XIV: u quark: bi, m, n 2 , di and p parameters corresponding to expression J6J when applied to uu channel for different WIMP 
masses. Mass independent parameters in JBJ for this channel are presented at the bottom of the table. 



Parameter 


WIMP mass interval (GeV) 


Fitting power law(s) 


bi 


50 < M < 8000 


2.96A/ UUMJb 


m 


50 < M < 8000 


2.91M" U3&i + 1.90M"- U1 ' X 


n 2 


50 < M < 8000 


0.0587M ui4b +0.848M^ u ' I4b 


di 


50 < M < 8000 


0.317M- uuduu + 0.403Af- U3Oi 


P 


50 < M < 8000 


4.74M u -« ay 



Table XV: Parameters corresponding to l|10[l for u quark. 61 and p parameters follow a simple power-law behavior in the whole studied 
WIMP mass interval. ri\, n 2 and d\ parameters follow a sum of two power laws in the whole mass interval. 



4- d quark 

For this channel there are five mass-dependent parameters: b\, ri\ ri2, c± and p. The mass independent parameters 
are a x = 5.20, b 2 = 5.10, d x = 0.410, c 2 = 0.0260, d 2 = 0.570 and q = 1.40 • 10 -4 . All parameters in this case are 
presented in Table HSU The mass range studied for this quark was from 50 to 5000 GeV. 

In this channel, no conclusion can be drawn about the existence of an asymptotic high mass limit in the spectral 
shape. The chosen parameters provide good fits from x = 2 ■ 10 -3 for M = 50 GeV whereas for the rest of masses the 
fits work very well till x — 5 • 10 -4 . 

The scaling of b\ with M is given by a simple power-law in the whole mass interval, whereas m, n 2 and c\ follow 
a sum of two power-law behavior in the whole studied mass. Finally, the p parameter presents a power-law behavior 
for M > 50 GeV. These results can be seen in Table IXVIII 



5. s quark 

For the s quark, there are just four mass dependent parameters b\, n 2 , d\ and p. The mass independent parameters 
for this particle in ([6]) are a\ = 4.83, ni = 2.03, b 2 = 6.50, c% — 0.335, c 2 — 0.0 (d 2 is irrelevant as for the u quark) and 
q = 2.40 • 10~ 4 . All these parameters are detailed in Table IXVlIII The studied mass range for this quark is between 
50 and 7000 GeV. 

As in the d quark case, no conclusion can be drawn about the existence of an asymptotic high mass limit in the 



WIMP mass (GeV) 


61 


m 


n 2 


Cl 


P 


50 


4.09 


2.69 


0.561 


0.327 


17.7 


100 


4.24 


2.47 


0.522 


0.293 


66.2 


200 


4.39 


2.34 


0.480 


0.258 


166 


500 


4.56 


2.28 


0.448 


0.228 


483 


1000 


4.75 


2.25 


0.426 


0.212 


1270 


2000 


4.91 


2.24 


0.409 


0.200 


3130 


5000 


5.11 


2.23 


0.394 


0.187 


10200 



\ai = 5.20 ; b 2 = 5.10 ; rii = 0.410 ; c 2 = 0.0260 ; d 2 = 0.570 ; q = 1.40 ■ 10~ 4 | 

Table XVI: d quark: b%, n\, n 2 , c\ and p parameters corresponding to expression (0 when applied to dd channel for different WIMP 
masses. Mass independent parameters in J6} for this channel are presented at the bottom of the table. 
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Parameter 


WIMr mass interval (GeV J 


Fitting power law(s) 


L 

01 


ou s ™ S ouuu 


Q QO »,fU.U4»b 


ni 


50 < M < 5000 


21 8M -u.uu3 + 2 .25M- u - uun;l 


112 


50 < M < 5000 


0.848A-f- ,Klu + 0.161Af uut>ra 


Cl 


50 < M < 5000 


0.722M~ U ' 2YU + 0.0544M UUS/4 


P 


50 < M < 5000 


0.168A/ l;iM 



Table XVII: Parameters corresponding to l|10|l for d quark. As can be seen, bi parameter follows a simple power law in the whole accessible 
mass interval. ri\ , n 2 and ci parameters follow a sum of two power laws in the whole accessible mass interval. Finally, p parameter follows 
a power law for M > 50 GeV. 



WIMP mass (GeV) 


bi 


n 2 


di 


P 


50 


4.78 


0.719 


0.367 


186 


100 


5.31 


0.669 


0.332 


409 


200 


5.43 


0.648 


0.315 


605 


500 


5.60 


0.612 


0.290 


1180 


1000 


5.73 


0.592 


0.276 


1980 


2000 


5.87 


0.575 


0.263 


3320 


5000 


6.07 


0.557 


0.249 


6500 


7000 


6.12 


0.548 


0.244 


7570 



ai = 4.83 ; m = 2.03 ; b 2 = 6.50 ; a = 0.335 ; c 2 = 0.0 ; q = 2.40 ■ 10~ 4 

Table XVIII: s quark: bi, n 2 , d\ and p parameters corresponding to expression JBJ when applied to ss channel for different WIMP 
masses. Mass independent parameters in JHJ for this channel are presented at the bottom of the table. 

spectral shape. The scaling with M of the parameters for this quark is a simple power law for b\ parameter for masses 
higher than 1000 GeV, the sum of two power laws for n,2 and d\ parameters in the whole studied WIMP mass, and 
two power laws for p parameter: one for masses smaller than 1000 GeV and another for masses higher than 1000 
GeV. These results are shown in Table IXIXl 



6. c quark 

As for the d quark, there are five mass dependent parameters. In this case b\, n\, ci, d\ and p which are presented in 
Table IXXl The mass independent parameters are a\ = 5.58, 62 = 7.9, = 0.686, C2 = 0.0 (therefore ^2 is irrelevant) 
and q = 9.00 • lO" 4 . 

Likewise the u quark, the studied mass range was from 50 to 8000 GeV and again no conclusion can be made 
about the existence of an asymptotic high mass limit for the spectral shape. Higher masses simulations would be thus 
required also in this case. 

The scaling of b\ and n\ with M shows a simple power-law behavior in the considered range. For ci and p, the 
single power-law evolution is only valid for masses above 200 GeV. Finally, di parameter follows a sum of two power 
laws in the studied mass range. These results are shown in Table IXXll 



Parameter 


WIMP mass interval (GeV) 


Fitting power law(s) 


bi 


M > 1000 


4.54M™ 


n,2 


50 < M < 7000 


3.68M- 1 " 1 + 0.744M~ U ' UJM 


di 


50 < M < 7000 


0.621M- Ut " 4 + 0.414M _UUO88 


P 


50 < M < 100 


4mM umi 




100 < M < 7000 


12.8M 0732 



Table XIX: Parameters corresponding to l|10[l for s quark. As can be seen, b± parameter follows a simple power-law behavior for masses 
higher than 1000 GeV. n 2 and d\ parameters follow the sum of two power laws for the whole studied WIMP mass interval. Finally, p 
parameter presents two power laws: one for masses smaller than 1000 GeV and another for masses higher than 1000 GeV. 
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WIMP mass (GeV) 


bi 


m 


Cl 


di 


P 


50 


5.93 


2.35 


0.239 


0.428 


210 


100 


5.48 


2.08 


0.283 


0.374 


379 


200 


4.98 


1.86 


0.330 


0.330 


673 


500 


4.50 


1.65 


0.378 


0.288 


1230 


1000 


4.00 


1.50 


0.406 


0.264 


2110 


2000 


3.70 


1.35 


0.432 


0.245 


4050 


5000 


3.27 


1.17 


0.470 


0.221 


8080 


8000 


3.08 


1.11 


0.494 


0.208 


12000 



|ai = 5.58 ; b 2 = 7.90 ; n 2 = 0.686 ; c 2 = 0.0 ; q = 9.00 ■ 10~ 4 | 

Table XX: c quark: b\, m, cj, d\ and p parameters corresponding to expression ((6J in the cc channel for different WIMP masses. Mass 
independent parameters in (0 for this channel are presented at the bottom of the table. 



Parameter 


WIMP mass interval (GeV) 


Fitting power law(s) 


bi 


50 < M < 8000 


9.90M- ui;!U 


m 


50 < M < 8000 


4.14M" u ' i4S 


Cl 


500 < M < 8000 


0.210M uuabi 


dx 


50 < M < 8000 


imM-" asA + 0.479M- UUU4 " 


P 


200 < M < 8000 


g 11M u.«i^ 



Table XXI: Parameters corresponding to l|10ll for c quark. It can be seen that the mass dependent parameters follow a power-law behavior 
for intermediate and high WIMP masses. In particular the d\ parameter follows a sum of two power-law behavior in he whole accessible 
WIMP mass range. 



7. b quark 

For the b quark, the required gamma rays spectra parametrization is the one given by expression (j6|). For this 
particle, the mass independent parameters are a\ = 10.0, b 2 = 11.0, c 2 = 0.0151, d 2 = 0.550, q — 2.60 • 10~ 4 . The 
mass dependent parameters are bi n\, n 2 , Ci, d\ and p. Their values are presented in Table IXXlTl 

The studied mass range is from 50 to 8000 GeV. Unlike previous particles for which the spectra did not change 
remarkably for very high masses, in the present case no conclusion can be drawn about the existence of an asymptotic 
high mass limit. 

Concerning the scaling behavior of the parameters for this quark, we observe that the behavior depends both on 
the WIMP mass and on the considered parameter. Thus b\ and n\ no longer scale with a single power-law for M 
higher than 100 GeV. For n 2 , two simple power laws can be seen, one from 50 to 1000 GeV (not included) and a 
second one from 1000 (included) to 8000 GeV. c\ shows also a power-law behavior but only up to 50 GeV. Finally, 
both d\ and p, scale with simple power laws from 500 GeV up. These results are summarized in Table IXXIIIl 

Some of these results are presented graphically in Figure 7, Appendix D for four WIMP masses: 50, 200, 1000 and 
5000 GeV . Also in this Appendix, mass dependent parameters bx, n%, n 2 , ci, d\ and p are plotted in Figure 8. 





WIMP mass (GeV) 


bi 


m 


n 2 


Cl 


di 


P 






50 


19.5 


6.48 


0.710 


0.365 


0.393 


57.8 






100 


17.1 


5.80 


0.695 


0.403 


0.360 


138 






200 


13.1 


5.01 


0.680 


0.415 


0.340 


281 






500 


8.76 


4.04 


0.660 


0.431 


0.319 


623 






1000 


6.00 


3.36 


0.647 


0.447 


0.305 


1030 






2000 


4.60 


2.85 


0.640 


0.460 


0.294 


1620 






5000 


3.00 


2.26 


0.634 


0.479 


0.280 


2670 






8000 


2.35 


2.00 


0.629 


0.490 


0.274 


3790 




ai = 


= 10.0 ; b 2 = 11.0 ; c 2 


= 0.0151 ; d 2 = 


0.550 


q = 2.60 • 10~ 4 



Table XXII: b quark: 6i, m, n 2 , ci, d± and p parameters corresponding to expression J6} in the bb channel for different WIMP masses. 
Mass independent parameters in JBJ for this channel are presented at the bottom of the table. 
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Parameter 


WIMr mass interval (tjeVJ 


Fitting power law(s) 


bi 


100 < M < 8000 


152Ai 


m 


100 < M < 8000 


18.7A/T"'™ 


n 2 


50 < Af < 1000 

i Ann ' Ti/r ^ onnn 

10UU < M < 8000 


n one a f—U Uiil.9 

0.805M 

n 7fi7 7\ /f - 0.0129 
0. (VIM 


Cl 


50 < A/ < 8000 


0.328Af UU44Y 


di 


50 < M < 600 
600 < M < 8000 


0.474M- u ' UMy +37.1Af" i ' 8Y 
0.449M- 0552 


P 


200 < M < 8000 


11.8M M41 



Table XXIII: Parameters corresponding to dlOI I for the b quark. All mass dependent parameters for the b quark follow simple power-law 
scalings at intermediate and high energies. Only n 2 and d\ parameters follow power-law behaviors at low WIMP masses. 

VII. NUMERICAL CODES 

All the calculations performed in this investigation are available at the website 
http:/ /teorica. fis.ucm.es/~PaginaWeb/downloads. html 

At this site, we provide the Mathematica [19] files that contain the fitting expressions ©, (|8j) and © for x lf3 diV 7 /d:r 
presented in this paper when applied for each studied channel. Let us remind that these parametrizations are valid 
in the corresponding WIMP masses intervals mentioned in the corresponding sections. Also in these files, the fitting 
formulae for mass dependent parameters in each channel are presented. 

VIII. CONCLUSIONS 

In this work, we have extensively studied the photon spectra coming from WIMP pair annihilation into SM particle 
particle-antiparticle pairs for all the phenomenologically relevant channels. The covered WIMP mass range has been 
optimized for each particular channel taking into account mass thresholds, statistics, and saturation of the Monte 
Carlo simulation. For instance, for light quarks it was from 50 GeV to 8000 GcV, for leptons it was from 25 GeV to 
50 TeV, for gauge bosons from 100 GeV to 1000 GeV and for t quark from 200 to 1000 GeV. All simulated spectra 
covered the whole accessible energy interval, from extremely low energetic photons till photons with one half of the 
available total center of mass energy. 

Once the spectra were obtained, analytical expressions were proposed to fit the simulation data. Three different 
fitting functions appeared to be valid depending on the studied channel: one for light quarks and leptons, another 
for gauge bosons and finally one for t quark very similar to the latter. Those expressions depended on either WIMP 
mass dependent or independent parameters. For WIMP mass independent parameters, their values did nevertheless 
depend on the considered annihilation channel whereas for WIMP mass dependent ones, their evolutions with WIMP 
mass were parametrized from the obtained values by continuous and smooth curves. 

In addition to a better understanding of the different channels for photon production from DM annihilation, the 
use of these fitting functions found in these analyses can save an important amount of computing time and resources: 
Monte Carlo simulations do not need to be repeated each time that a particular photon spectrum needs to be known 
for a given channel and center of mass energy. This fact is particularly important for high energy photons, whose 
production rate is very suppressed and would require large computation times and to store big amounts of data. Our 
research was thus able to present very good statistics for those energies. 

By having used extensive PYTHIA Monte Carlo simulation we have been able to obtain relatively simple 
parametrizations of these spectra and fit the corresponding parameters. As our analysis is model independent, it 
could be useful, both for theoreticians and experimentalists, interested in the indirect DM detection through gamma 
rays. Given some theoretical model, and the corresponding velocity averaged annihilation cross sections for the dif- 
ferent channels, our formulae make it possible to obtain the expected photon spectrum for each particular theoretical 
model in a relatively simple way. In this sense, further work is in progress to extend our analysis to other stable 
particles like positron or neutrinos but these results will be presented elsewhere. 
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IX. APPENDICES 

In this section we present simulations for some of the studied channels: W + W~, ti, t + t~ and bb. For these 
channels four simulated spectra are presented together with the proposed fit formulae. For each channel, evolution 
with WIMP mass of mass dependent parameters have been plotted. The final appendix E shows the running with 
the WIMP mass of the total number of photons per WIMP pair annihilation. 



A. Plots for W gauge boson 




(a) Photon spectrum for M = 100 GeV for W + W channel. (b) Photon spectrum for M = 200 GeV for W + W channel. 




(c) Photon spectrum for M = 350 GeV for W + W channel. (d) Photon spectrum for M = 1000 GeV for W + W channel. 



Figure 1: Photon spectra for four different WIMP masses (50, 200, 1000 and 5000 GeV) in the W + W~ channel. Red dotted 
points are PHYTIA simulations and solid lines correspond to the proposed fitting functions. 
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Figure 2: Mass dependence of bi, ci, di, p and j parameters for W + W channel. Crossed points are parameters values found 
after the fitting process for each WIMP mass and solid lines correspond to the proposed fitting functions. 
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Figure 3: Photon spectra for four different WIMP masses (200, 250, 500 and 1000 GeV) in the tt annihilation channel. Red 
dotted points are PHYTIA simulations and solid lines correspond to the proposed fitting functions. 
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(b) ni parameter of expression J6j for tt channel. 



Fitting power law 
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(c) C2 parameter of expression J6jl for tt channel. 



Fitting power law + Smooth spline 
p values tor top-anti top channel 
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WIMP mass (GeV) 

(d) p parameter of expression JHJ for tt channel. 



Fitting power law 
s tor top-anti top channel 



400 500 



WIMP mass (GeV) 

(e) q parameter of expression JH} for tt channel. 



Fitting power law 
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(f) / parameter of expression (6]) for tt channel. 



Figure 4: Mass dependence of fri, ni, C2, p, q and I parameters for tt annihilation channel. Crossed points are parameters values 
found after the fitting process for each WIMP mass and solid lines correspond to the proposed fitting functions. 
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C. Plots for r lepton 



Total fit 
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Total fit ■ 
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(a) Photon spectrum for M = 25 GeV for r + r channel. (b) Photon spectrum for M = 100 GeV for r + r channel. 
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(c) Photon spectrum for M = 1000 GeV for r+r channel. (d) Photon spectrum for M = 5 • 10 4 GeV for r+r channel. 

Figure 5: Photon spectra for four different WIMP masses (25, 100, 1000 and 5 ■ 10 4 GeV) in the r + r" annihilation channel. 
Red dotted points are PHYTIA simulations and solid lines correspond to the proposed fitting functions. 



Fitting power laws 
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Fitting power laws - 
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WIMP mass (GeV) 

(a) ni parameter of expression {6} for r + r~ channel. 



WIMP mass (GeV) 

(b) p parameter of expression {6jl for t + t — channel. 



Figure 6: Mass dependence of ni and p parameters for r + r annihilation channel. Crossed points are parameters values found 
after the fitting process for each WIMP mass and solid lines correspond to the proposed fitting functions. 
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D. Plots for b quark 
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(a) Photon spectrum for M = 50 GeV for bb channel. 



E y /M W!MP 



(b) Photon spectrum for M = 200 GeV for bb channel. 




(c) Photon spectrum for M = 1000 GeV for bb channel. 



E y /M W!MP 

(d) Photon spectrum for M = 5000 GeV for bb channel. 



Figure 7: Photon spectra for four different WIMP masses (50, 200, 1000 and 5000 GeV) in the bb annihilation channel. Red 
dotted points are PHYTIA simulations and solid lines correspond to the proposed fitting functions. 
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Figure 8: Mass dependence of bi, m, n2, ci, d± and p parameters for the bb annihilation channel. Crossed points are parameters 
values found after the fitting process for each WIMP mass and solid lines correspond to the proposed fitting functions. 
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E. Photon number per WIMPs annihilation 




Figure 9: Total photon number per WIMP annihilation: On the left (a) the leptonic channels are presented whereas on the 
right (b) both gauge bosons and quarks results are plotted. 



